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New measurements are re-ported on the solubility of germanium in liquid 
gallium, thallium, tin, arsenic, bismuth, cadmium and zinc, and the solu- 
bility of silicon in liquid indium, tin, lead, antimony, bismuth and zinc. 
The measurements of other workers are reviewed, including those of the 
solubility of germanium and silicon in liquid copper, silver, gold and alumi- 
num; of germanium in liquid indium, lead and antimony; and of silicon 
in liquid arsenic and nickel. All but two of the liquidus curves can be de- 
scribed within experimental error by a two-constant equation. The form of 
this equation suggests that the liquid solutions exhibit certain simple ther- 
modynamic properties, and some evidence is cited indicating that the con- 
stants of this equation can be used to estimate the excess free energy of the 
solutions. Figures for the complete liquidus curves of these binary systems 
(T-x and log x-l/T) have not been included in this paper, but sets of these 
figures can be supplied upon request. 

I. INTRODUCTION 

Germanium and silicon have become very important elements in the 
last ten years. Great strides have been made in the exploitation of the 
electrical properties of these semiconductor elements, and many use- 
ful solid state electronic devices that employ these two elements are 
manufactured in increasing numbers every year. The electrical proper- 
ties of pure germanium and silicon are of academic interest only, how- 
ever. It is the modification of these properties by small amounts of added 
impurities that has made these elements technologically important. 

A knowledge of the chemical reactivity of electrically active impuri- 
ties with germanium or silicon is necessary in order to be able to control 
the impurity concentrations. Much of the chemistry of such reactions 
can be described in a simple way by the use of phase diagrams. These 
diagrams give the chemical compositions of phases at equilibrium at 
various temperatures. Of particular interest are the liquidus and solidus 
curves of germanium and silicon. 
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The liquidus curve gives the solubility of germanium or silicon in 
another element as a function of temperature. Each point on the curve 
corresponds to the temperature and composition at which a liquid phase 
is saturated with germanium or silicon respectively. The solidus curve 
gives the concentration of the other component which will be dissolved 
in germanium or silicon at various temperatures when the solid phase 
is in equilibrium with a saturated liquid phase. 

In this paper, the liquidus curves of germanium and silicon will be 
discussed; the solidus curves are discussed in an accompanying paper. 1 

The liquidus curves formed by the dissolution of germanium and sili- 
con in elements of Groups III, IV and V of the periodic system were 
reviewed a number of years ago. 2 It was found that these liquidus curves 
could be reasonably well approximated by a one-constant equation that 
had a form consistent with the thermodynamic properties expected of a 
"regular" liquid solution in equilibrium with a pure solid phase. Addi- 
tional solubility measurements now show that the liquidus curves can- 
not be satisfactorily represented by such an equation. A simple two- 
constant equation has been found that describes the curves within the 
present limits of experimental error. This equation suggests that simple 
departures from ideal solution entropy occur, as well as a simple heat-of- 
mixing effect. 

Our objectives have been: (a) to construct the best liquidus curves of 
germanium and silicon from measurements reported in the literature, 
supplemented by our own measurements and (b) to obtain the activity 
coefficients of germanium and silicon along the liquidus curves. This in- 
formation will be used to estimate the thermodynamic properties of the 
liquid binary alloys. 

One method of evaluating liquidus curves is to plot the solubility 
measurements on a temperature-composition (atom fraction) scale and 
draw the best curve through the points from the eutectic to the melting 
point. For many of the liquidus curves of interest here this is neither 
practical nor desirable. A number of the liquidus curves cover a very long 
temperature range (as much as 1380 degrees in the case of the silicon- 
gallium system) and long composition range (from an atom fraction of 
about 10~ 10 to 1 in the same system). There is only a limited amount of 
experimental data available for these systems and rather long extrapola- 
tions of the data, to both higher and lower temperatures and composi- 
tions, are required if an estimate of the complete curve is to be made. 

It is possible, however, to make such extrapolations with greater con- 
fidence by taking advantage of all other available information that in- 
fluences the position of a liquidus curve. The first and most important 
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bit of additional information is that the measured liquidus curves appear 
to represent equilibrium conditions between liquid and solid phases. 
Studies of solubility as a function of time lead to this conclusion, as do 
the agreement between measurements made by different methods. The 
disagreement that exists between the measurements of solubility made 
by different workers frequently arises, Ave believe, from experimental 
error, part of which may be attributable to the fact that thermodynamic 
equilibrium was not attained. However, in general, this problem is 
avoidable. 

With the knowledge that thermodynamic equilibrium exists between 
the phases, additional information becomes pertinent: (a) germanium 
and silicon exhibit only one crystalline modification; (b) the solid solu- 
bility of the other component in germanium and silicon is small; (c) the 
melting points and heats of fusion of germanium and silicon are now 
better known; (d) the liquidus curves can be expected to exhibit certain 
simple properties in the region of the melting point of germanium and 
silicon; (e) the liquidus curves that extend to low temperatures and 
represent liquids that become very dilute in germanium or silicon can 
also be expected to exhibit certain simple properties. 

The implications of this additional information are as follows: (a) the 
fact that germanium and silicon exhibit only one crystalline modifica- 
tion means that the liquidus curves will be smooth curves with no abrupt 
changes of slope; (b) the fact that the solid solubility of the second 
component is small means that the thermodynamic properties of the 
solid phase are essentially those of pure germanium or silicon; (c) the 
more accurate knowledge we now have of the melting points of ger- 
manium and silicon can be used to construct better curves, since curves 
now in the literature use melting points varying from 936 to 958°C for 
germanium and from 1400 to 1430°C for silicon; (d) a knowledge of the 
heats of fusion of germanium and silicon and the fact that the solid solu- 
bilities are small means that the liquidus curves can be expected to come 
into the melting point of germanium or silicon with a certain known 
limiting slope; (e) when the amount of solute in the liquid phase is small 
at low temperatures, it can be expected that the logarithm of the atom 
fraction of the solute will be well approximated by a straight line when 
it is plotted as a function of the reciprocal of the absolute temperature. 

We have used a method of evaluating the liquidus curve measure- 
ments that takes advantage of all the above information related to these 
curves. In addition, this method of evaluation has led to the discovery 
of an apparent regularity that leads to an interesting suggestion about 
the thermodynamic properties of the liquid alloys. Specifically, we have 
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calculated the activity coefficients of germanium and silicon along a 
liquidus curve from the solubility measurements, and calculated the 
parameter a = RT In 71/ (1 — x) , where R is the gas constant, T the 
absolute temperature, 71 the activity coefficient of germanium or silicon 
referred to the pure supercooled liquid state, and x the atom fraction of 
germanium or silicon in the saturated liquid. It has been found that a 
is a linear function of T. This leads to a two-constant liquidus curve 
equation that has all the desired features. Previously 2 it had been found 
that the liquidus curves of a number of these systems could be approxi- 
mated by a one-constant equation, but more accurate solubility meas- 
urements now show that two constants are needed. 

The liquidus curves formed by 14 different elements with germanium 
and silicon, for a total of 27 different curves, have been evaluated. The 
parameter a is found to be a linear function of temperature over a sig- 
nificant temperature range within experimental error for 17 of these 
binary systems (germanium-indium, -gallium, -aluminum, -lead, -tin, 
-bismuth, -cadmium, -silver, -copper, -gold; silicon-indium, -gallium, 
-itin, -bismuth, -antimony, -silver, -zinc). The solubility data for seven 
other systems (germanium -arsenic, -antimony; silicon-aluminum, -lead, 
-arsenic, -copper, -gold) are more limited or exhibit considerable scatter. 
An estimate of the best liquidus curve representing the available meas- 
urements is obtained by assuming that a is a linear function of T. The 
complete germanium-thallium liquidus curve has been estimated from 
very limited solubility measurements, and the silicon-thallium curve 
has been estimated, although no solubility measurements have been 
made. Two systems, germanium-zinc and silicon-nickel, have liquidus 
curves that appear to be qualitatively the same as the other curves 
when plotted as T versus x, but it is found that a is not a linear func- 
tion of temperature. New solubility measurements are reported for 14 
binary systems (germanium-thallium, -gallium, -tin, -bismuth, -cad- 
mium, -arsenic, -zinc; silicon-tin, -lead, -bismuth, -antimony, -zinc and 
the previously unreported work of Hassion 3 on germanium-indium and 
germanium-lead ) . 

After a description of our experimental procedure we will illustrate 
this method of evaluation of the solubility measurements by discussing 
the germanium -gallium system in some detail. Each of the other sys- 
terms will then be considered briefly. 

II. EXPERIMENTAL PROCEDURE 

The solubility measurements were made by a method similar to that 
used by Kleppa and Weil. An excess of germanium or silicon was sealed 
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into an evacuated silica tube with a known weight of solvent. The ger- 
manium and silicon were in the form of single-crystal ingots of high 
purity. The temperature of the furnace containing the quartz tube was 
raised slowly to insure against overshooting, and then held constant for 
various periods of time before the furnace was tipped. A constriction in 
the tube permitted separation of the saturated melt from the excess 
germanium or silicon by tipping of the furnace. The equilibration time 
varied from one-half to one hour for germanium and up to 19 hours for 
silicon. During equilibration the temperature was held constant to within 
one degree centigrade, by manual control for the shorter times and by a 
Leeds and Northrup Speedomax controller for the longer heating times. 
The emf's of the calibrated platinum and platinum-10 per cent rhodium 
thermocouples used for temperature measurement were determined with 
a Leeds and Northrup portable precision potentiometer. Temperatures 
are considered to be accurate to within ±2 degrees. 

The compositions of the saturated melts were obtained from the loss 
in weight of the germanium or silicon ingots, except for the germanium- 
arsenic system, in which case the poured-off melts were chemically 
analyzed. In order to obtain the loss in weight of the single-crystal in- 
gots, the small amount of melt adhering to the crystal was dissolved in 
a solvent that would not attack the pure germanium.* 

III. TREATMENT OF DATA 

3.1 Germanium-Gallium 

The solubility of germanium in liquid gallium has been measured by 
Klemm et al., 5 Keck and Broder, 6 Greiner 7 and de Roche. 8 These meas- 
urements are shown in Fig. 1, along with our own measurements. The 
measurements of Klemm et al. are in substantial disagreement with the 
others. The disagreement between the measurements of de Roche and 
the others at the lower temperatures is quite evident in Fig. 2, where 
log x has been plotted as a function of the reciprocal of the absolute 
temperature, 1/T. The measurements of Klemm et al. are not included 
in this figure. 

An interesting feature of these measurements is that they fall close to 
the ideal liquidus curve. Included in Fig. 2 are two dashed curves, each 
an ideal liquidus curve corresponding to different assumptions about 
the heat capacity of supercooled liquid germanium. One of the curves 

* In general, hot concentrated HC1 was used. A mixture of 30 per cent H2O2 , 
glacial HAc and H«0 in the volume ration 2:2:5 was used for lead on germanium. 3 
Aqua regia was used for antimony on silicon. 
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Fig. 1. — Germanium-gallium liquidus curves. 
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Fig. 2 — Germanium-gallium liquidus curves, log x vs. 1/T. 
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results from the assumption that the heat capacity of liquid germanium 
is the same as solid germanium (AC P = 0); the other results from the 
assumption that the heat capacity of liquid germanium is constant 
(AC P ^ 0). If the liquid alloys are ideal solutions of atoms, if these 
liquids are at equilibrium with pure solid germanium and if the heat 
capacity of pure supercooled liquid germanium is the same as pure solid 
germanium, the liquidus curve will be given by the following equation 
(Ref. 9, pp. 357-358): 



, AH/ 

In x = 



h-f} (1) 



R \Ti 

If, however, the heat capacity of pure supercooled liquid germanium is 
constant and the heat capacity of pure solid germanium has the form 
a + bT, the ideal liquidus curve will be given by 



-=¥fe-»+M<— > 



Ti° - T , , T 
~^r- + a]n- 

+ *b[CTi ) 1 - T 2 ] 



(2) 



The dashed curve of Fig. 2 designated AC P = corresponds to ( 1 ) 
with AH X F , the latent heat of fusion of germanium, taken as 8100 calories 
per gram atom 10 and Ti°, the melting point of germanium, 937°C. 7,11 
The dashed curve of Fig. 2 designated AC P ^ is obtained from (2), 
using the specific heat equation listed by Kubaschewski and Evans " for 
solid germanium and assuming that the specific heat of liquid germanium 
is constant and equal to that of solid germanium at its melting point. 

At temperatures near the melting point of germanium, de Roche re- 
ports that the lowering of the freezing point leads to a heat of fusion of 
8100 ± 200 calories per gram atom. This is in agreement with Greiner 
and Breidt. 10 The agreement between the heat of fusion obtained by 
freezing point lowering measurements and the heat obtained by direct 
measurement is evidence that germanium and gallium form atomically 
dispersed liquid solutions.* If other elements form atomically dispersed 
binary liquid phases with germanium and also exhibit negligible solid 
solubility, their liquidus curves will all come into the melting point of 
germanium with the same limiting slope. In other words, (1) will be a 
good approximation to all the liquidus curves at temperatures near the 
melting point of germanium. 



* The possibility cannot be ruled out that germanium and gallium form poly- 
meric species (Cle)* and (Ga)i , but this is believed to be unlikely. 
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At low temperatures, where the atom fraction of germanium in the 
saturated liquid phase is small, we expect that the form of (1) will be 
correct but that, in general, the slope and intercept (on a log x-l/T 
plot) will be different for each element. 

Let it now be assumed that ( 1 ) is the correct equation for the ideal 
liquidus curve (this corresponds to saying that the heat of fusion of 
germanium is not a function of temperature — i.e., the heat capacity of 
solid and supercooled liquid germanium are the same — and that the solid 
solubility of the component forming the ideal liquid solutions with ger- 
manium is zero). We can measure the departures of any experimental 
liquidus curve from the ideal liquidus curve in terms of the activity 
coefficient, 71 , defined by the expression (Ref. 9, pp. 357-358) : 



, AH 

In ZYi = — ^ 



l -Q.-i)- w 



Consequently, a value of 71 can be calculated for every solubility meas- 
urement. 

The activity coefficient defined in this way has the pure (supercooled) 
liquid germanium as the reference state. Thus, when x is unity, 71 is 
unity. Since x will be unity at the melting point of germanium, ( 1 ) be- 
comes a good approximation to (3) when the temperature is near the 
melting point of germanium. 

The relationships between 71 and other thermodynamic parameters 
are given below : 

RT In 71 = Aft', (4) 

AFi = Mix - TASi. (5) 

The terms Af\ e and ASi are the relative partial molar excess free energy 
and entropy, respectively, and A#i is the relative partial molar enthalpy. 
From (4) and (5), and using the fact that the heat of fusion is equal to 
Ti°ASi F , (3) may be rewritten in the following form: 

AHi" + A#i . ASi F + ASi ,„, 

lnx= w + g . (6) 

When the liquidus curve extends to low temperatures where x is 
small, AHi and ASi will no longer depend significantly on composition. 
Since the temperature dependences of AHi and ASi are expected to be 
small, (6) will give a linear relationship between In x and 1/T. In the 
germanium-gallium system, the values of AHi and ASi are close to 
zero. 
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Equation (6) can also be used as a general expression for the liquidus 
curves, where, in general, AH and ASi are functions of composition and 
temperature. Since we expect the temperature dependence of AR X and 
A*Si" to be small, the variation in AH l and ASi will arise primarily from 
their composition dependence. The composition dependence of Ai? x and 
A*Si" can each be represented conveniently by a power series in (1 — x) : 

AH, = £a„(l - x) n , (7) 

A&' = EMI -x)\ (8) 

Since Ai?i and ASi are zero when x = 1, a and 6 are both zero. If 
long-range forces between the atoms are not present, a, and &i will be 
zero. 13 Consequently, a first approximation to the composition de- 
pendences of AHi and ASi" can be expected to be 

A/?! = a 2 (l - x? (9) 

and 

A& e = 6 2 (1 - x)\ (10) 

Upon substitution of (9) and (10) in (6) and rearranging, we obtain 

TAS" - AH F - RT In x , ^ n 1 , 

a = r- r^ = a 2 — Oil . (11; 

(i - x)- 

A value of the parameter a can be calculated from every solubility meas- 
urement. If A/7 and AS" are independent of temperature and their com- 
position dependences are given by (9) and (10), it follows that a will 
be a linear function of T. 

We have evaluated a for a number of the liquidus curves of germanium 
and silicon and found that the experimental data can be satisfactorily 
approximated with a as a linear function of T. As will be discussed below, 
it does not follow that (9) and (10) must represent the thermodynamic 
properties of the liquid phases in order for a to be a linear function of T. 
For any given set of experimental points there are an infinite number of 
equations for Ai? x and A^i" that would provide a fit to an a-T plot within 
experimental error. The discussion of the preceding paragraph describes 
the simplest thermodynamic properties the liquid alloys could have 
which would lead to the liquidus curve shapes we observe. A previous 
study 2 of a number of these liquidus curves suggested that a was inde- 
pendent of temperature along the liquidus curve. This property of a 
would be expected if the liquid alloys were regular solutions.* 

* The term "regular solutions" was coined b} r Hildebrand to denote solutions 
for which "thermal agitation is sufficient to give practically complete random- 
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The solubility measurements of germanium and gallium evaluated in 
this way are shown in Fig. 3. The measurements of Greiner, Keck and 
Broder 6 and our measurements can be represented by a line with a = 
— 150 and 6 = 0. The measurements of de Roche 8 are included, except 
those at low temperatures and those near the melting point of germa- 
nium. The parameter a is quite sensitive to experimental error near the 
melting point, and, between about 1100 and 1211°K, a = —150 is a 
satisfactory representation of all the measurements when a T-x plot is 
used. 

The solubility measurements we have made of germanium in gallium 
are given in Table I ; the measurements of Keck and Broder, scaled from 
their published figure, and of Greiner are also given. These measurements 
have been used to obtain a and b. 

3.2 Germanium-Aluminum, Germanium-Indium 

The solubility of germanium in aluminum has been measured by 
Stohr and Klemm. 16 Their results, scaled from a figure, are given in Table 
I.f The parameter a has been plotted as a function of temperature in 
Fig. 3. A line, a = -5360, 6 = -3.16 has been drawn through the 
points. Some curvature is suggested, but it is believed that the straight 
line fits the data within experimental error. 

The solubility of germanium in indium has been measured by Keck 
and Broder, 6 whose data, scaled from a figure, are given in Table I, 
and by Hassion, 3 whose measurements are also given in the table. The 
measurements of Klemm et al. 5 are in sufficient disagreement with these 
measurements to justify not including them in Table I and Fig. 3. A 
line, a = 1570, b = 0.56, has been drawn through the measurements 
of Keck and Broder and Hassion, plotted as a versus T in Fig. 3. 

3.3 Germanium-Tin, Germanium-Lead 

The solubility of germanium in tin has been measured by Stohr and 
Klemm. 16 However, we have also measured the solubility of germanium 

ness. 14 This definition would, in general, include those solutions for which A»Si c is 
zero and for which (9) gives the partial molar heats of solution, but would not be 
restricted to them. It is useful to classify solutions empirically in terms of the 
mathematical functions used to describe their thermodynamic properties. The 
terms "regular solutions" and "strictly regular solutions" are sometimes used 
(Ref. 9, p. 246 and Ref. 15, p. 85) to describe solutions for which a 2 is constant 
and b« is zero in (9) and (10). 

t The parameter a is very sensitive to error in T and x near the melting point. 
Consequently, solubility measurements in this range have not been weighted very 
heavily and frequently have not been included in the table of data. 
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in tin and are in significant disagreement. Our measurements are given 
in Table I and a plot of a versus T is shown in Fig. 4. A line, a = 1680, 
b — 1.08, has been drawn through the points. The measurements of 
Stohr and Klemm when plotted as a versus T fall on a curve which 
crosses our data at about 800°K and is of opposite slope. Thus, around 
800°K, our measurements agree with those of Stohr and Klemm, but at 
higher temperatures we measure somewhat higher solubilities and at 
lower temperatures somewhat lower solubilities. 
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Table I — Experimental Compositions (Atom Fraction 

of Germanium) and Temperatures of Germanium 

Liquidus Curves 



Element 


Source 


Temperature, 
"K 


X 


Aluminum 


Stohr & Klemm (Ref. 16) 


762 


0.324 






849 


0.410 






926 


0.501 






1014 


0.614 






1081 


0.718 






1172 


0.904 


Gallium 


This work 


570 


0.0268 






653 


0.0655 






749 


0.142 






749 


0.144 






843 


0.243 






933 


0.382 




Keck & Broder (Ref. 6) 


470 


0.0042 






473 


0.0052 






470 


0.0055 






570 


0.026 






670 


0.071 




Greiner & Breidt (Ref. 7) 


1018 


0.541 






1134 


0.796 






1171 


0.894 






1191 


0.944 


Indium 


Keck & Broder (Ref. 6) 


575 


0.0061 






575 


0.0083 






683 


0.033 






769 


0.076 




Hassion (Ref. 3) 


735 


0.064 






735 


0.068 






756 


0.073 






830 


0.138 






830 


0.133 






862 


0.160 






907 


0.254 






921 


0.249 


Thallium 


This work 


1023 


0.131 






1073 


0.260 






1123 


0.656 




Klemm et al. (Ref. 5) 


1052 


0.238 






1106 


0.377 






1122 


0.492 






1143 


0.635 






1165 


0.776 


Tin 


This work 


601 


0.014 






658 


0.0305 






758 


0.077 






810 


0.125 






860 


0.181 






908 


0.257 






960 


0.362 






999 


0.459 
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Table I — Continued 



Element 


Source 


Temperature, 
°K 


X 


Lead 


Hassion (Ref. 3) 


901 


0.021 






953 


0.040 






984 


0.052 






1012 


0.066 






1050 


0.105 






1058 


0.112 






1062 


0.126 






1108 


0.270 






1123 


0.37 






1137 


0.50 






1146 


0.63 


Antimony 


Ruttewit & Masing (Ref. 18) 


872 


0.20 






907 


0.25 






977 


0.40 






1043 


0.55 






1058 


0.60 






1100 


0.70 


Arsenic 


This work 


1023 


0.595 






1073 


0.693 






1123 


0.783 






1148 


0.841 




Stohr & Klemm (Ref. 16) 


1023 


0.605 






1102 


0.738 






1128 


0.798 


Bismuth 


This work 


873 


0.0273 






923 


0.0440 






973 


0.0656 






973 


0.0728 






1023 


0.118 






1073 


0.201 


Copper 


Reynolds & Hume-Rothery (Ref. 21) 


917 

948 

982 


0.365 
0.407 
0.460 










1010 


0.507 


r 




1048 


0.583 






1090 


0.674 






1154 


0.833 


Silver 


Maucher (Ref. 20) 


925 


0.241 






929 


0.246 






951 


0.271 






1028 


0.389 






J 104 


0.598 






1131 


0.690 






1151 


0.774 




Briggs et al. (Ref. 22) 


923 


0.260 






948 


0.270 






953 


0.286 






978 


0.308 






988 


0.328 






1013 


0.357 






1026 


0.398 
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Table I — Concluded 



Element 


Source 


Temperature, 
°K 


X 






1070 


0.497 






1103 


0.618 






1128 


0.708 






1149 


0.803 






1186 


0.923 


Gold 


Jaffee et al. (Ref. 23) 


629 


0.270 






919 


0.537 






1070 


0.731 






1164 


0.891 


Zinc 


This work 


723 


0.0914 






773 


0.148 






823 


0.213 






873 


0.288 






928 


0.368 






973 


0.440 






1023 


0.530 






1073 


0.636 


Cadmium 


This work 


669 


0.0076 






669 


0.0076 






722 


0.0144 






776 


0.0284 






776 


0.0282 






776 


0.0292 






821 


0.0469 






877 


0.0717 






927 


0.131 






978 


0.205 






1015 


0.306 






1015 


0.298 






1015 


0.279 



The solubility of germanium in lead has been measured by Briggs and 
Benedict, 17 Ruttewit and Masing 18 and Hassion. 3 We have accepted the 
measurements of Hassion given in Table I, and have plotted these data 
as a versus T, shown in Fig. 4. A line, a = 8780, b = 4.08, has been 
drawn to represent these measurements. The measurements of Briggs 
and Benedict are in fail- agreement with Hassion's measurements, but 
the low-temperature solubilities reported by Ruttewit and Masing are 
much lower than the values obtained from an extrapolation of Hassion's 
measurements using the a-T plot. 

3.4 Germanium- Arsenic, Germanium- Antimony, Germanium-Bismuth 

The measurements of Stohr and Klemm 16 of the solubility of ger- 
manium in arsenic have been scaled from their published figure and are 
listed in Table I along with our measurements. The a-T plot is shown 
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in Fig. 5 and a line, a = —5600, b = —4.16, has been drawn to repre- 
sent the data. 

The solubility of germanium in antimony has been measured by Rut- 
tewit and Masing 18 and Stohr and Klemm. 16 The measurements of Rut- 
tewit and Masing, scaled from a figure, are given in Table I and the 
a-T plot resulting from these measurements is shown in Fig. 5. The line, 
a = 2640, b = 1.98, has been drawn. Stohr and Klemm found smaller 
solubilities at all temperatures than those found by Ruttewit and 
Masing. A plot of a versus T with their data leads to a curve lying 500 
to 1,000 calories higher than the line drawn to represent the measure- 
ments of Ruttewit and Masing. We consider the antimony-germanium 
liquidus curve to be in some doubt. 

The solubility of germanium in bismuth has also been measured by 
Ruttewil and Masing 18 and by Stohr and Klemm. 16 We have also meas- 
ured the solubility of germanium in bismuth, and our results are given 
in Table I and plotted in Fig. 5. The line, a = 5505, b = 1.49, has been 
used to represent the data. Neither the measurements of Ruttewit and 
Masing nor those of Stohr and Klemm can be represented as a line on 
a plot of a versus T. Their measurements are in only fair agreement 
with those we have made. 

3.5 Germanium-Copper, Germanium- Silver , Germanium-Gold 

The solubility of germanium in copper has been measured by Schwarz 
and Elstner, 19 M aucher 20 and Reynolds and Hume-Rot hery. 21 The meas- 
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Fig, 5 — a vs. T for germanium-arsenic, germanium-antimony and germanium- 
bismuth. 

urements of Reynolds and Hume-Rothery are given in Table I and the 
corresponding a-T plot is shown in Fig. 6. The measurements of the other 
workers scatter rather badly. They are in agreement with Reynolds and 
Hume-Rothery at low temperatures but show lower solubilities at higher 
temperatures. The line, a = —7360, b = —7.67, has been drawn in 
Fig. 6. 

The solubility of germanium in silver, measured by Maucher, and 
by Briggs, McDuffie and Willisford, 22 is given in Table I; the a-T plot 
of these data is shown in Fig. 6. The line, a — —5500, b = —7.13, has 
been drawn. 

The solubility of germanium in gold has been measured by Jaffee, 
Smith and Gonser. 23 Their measurements are given in Table I and the 
a-T plot in Fig. 6. The line, a = —4865, b = — 1.02, has been drawn. 



3.6 Germanium-Zinc, Germanium-Cadmium 

Gebhardt 24 and Kleppa and Thalmayer 25 have measured the solubility 
of germanium in zinc. The measurements of the latter are in fairly good 
agreement with our measurements, which are recorded in Table I. It 
can be seen in Fig. 7 that a versus T cannot be represented by a straight 
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line. The measurements of Gebhardt and of Kleppa and Thalmayer are 
included in the plot. 

We have measured the solubility of germanium in cadmium, and the 
results are given in Table I. The points of the a-T plot appear to be 
represented within experimental error by the line, a = 41 10, 6 = 1.75. 
Spengler 26 has reported that the germanium-cadmium system is simple 
eutectic and has mentioned three solubility measurements. These three 
points are at much higher germanium concentrations than those we have 
found. 

3.7 Silicon-A luminum, Silicon-Gallium, Silicon-Indium 

We have evaluated the measurements of the .solubility of silicon in a 
number of elements of interest in a manner completely analogous to that 
used for germanium. We have used a heat of fusion of 12,100 calories 
per gram atom. 27 The melting point of silicon has been reported as 
141 5°C by Gayler 28 and 1408°C by Pell. 29 We have used a value of 
1410°C. 

The solubility of silicon in aluminum has been measured by Fraenkel, 
Roberts, 31 Gwyer and Phillips, 32 Broniewski and Smailowski 33 and Craig- 
head, Cawthorne and Jaffee. The measurements of Fraenkel and 
Broniewski and Smailowski were scaled from the figure published by 
Hansen. 35 These measurements are listed in Table II and a versus T is 
plotted in Fig. 8.* The line, a = —4140, b = —1.22, has been drawn. 

The silicon-gallium liquidus curve has been studied by Keck and 
Broder and Klemm et al. The measurements tabulated in Table II were 
scaled from figures. The points plotted in the a-T plot of Fig. 8 have 
been represented by the line, a = 3250, b = 0.83. 

The solubility of silicon in indium has been measured by Keck and 
Broder 6 and by Klemm et al. We have also studied this system. Our 
solubility measurements are given in Table II and the a-T plot is shown 
in Fig. 8, where the line, a = 11,450, b = 3.37, has been draAvn. The meas- 
urements of Klemm et al., are in essential agreement with our measure- 
ments, but by themselves give no indication of a linear relation between 
a and T. The three measurements of Keck and Broder are all at higher 
solubilities than we have measured, their lowest temperature measure- 
ment being in greatest disagreement with the liquidus curve correspond- 
ing to the line of Fig. 8. 

* The tabulation in Table II does not include the measurements made near the 
melting point of silicon. The measurements of Gwyer and Phillips near the eutec- 
tic, which are in agreement with the measurements of Craighead et al., are not 
plotted in Fig. 8. 



Table II — Experimental Compositions (Atom Fraction 

of Silicon) and Temperatures of Silicon 

Liquidus Curves 



Element 


Source 


Temperature, 
°K 


X 


Aluminum 


Fraenkel (Ref. 30) 


957 


0.194 






1232 


0.391 






1439 


0.597 






1503 


0.717 




Roberts (Ref. 31) 


969 


0.179 






1116 


0.319 






1232 


0.391 






1307 


0.482 






1503 


0.698 






1526 


0.731 






1607 


0.850 




Broniewski & Smailowski (Ref. 33) 


1131 
1333 


0.337 
0.510 




Gwyer & Phillips (Ref. 32) 


850 


o!ll3 






852 


0.116 






852 


0.119 






867 


0.124 






875 


0.134 






889 


0.142 






950 


0.183 




Craighead et al. (Ref. 34) 


851 


0.118 






854 


0.123 






871 


0.132 


Gallium 


Keck & Broder (Ref. 6) 


548 


0. 00001 6 






769 


0.00254 






1023 


0.034 






1023 


0.050 






1170 


0.100 






1273 


0.200 






1423 


0.400 


Indium 


This work 


1173 


o.ooss 






1273 


0.0190 






1273 


0.0205 






1373 


0.0425 






1373 


0.0494 






1373 


0.0400 






1473 


0.0955 






1473 


0.101 






1473 


0.101 


Tin 


This work 


1025 


0.00391 






1073 


0.00607 






1075 


0.00792 






1075 


0.00760 






1099 


0.00754 






1099 


0.00778 






1150 


0.0112 






1150 


0.0113 






1173 


0.0143 






1173 


0.0136 






1202 


0.0166 






1223 


0.0258 



1S7 
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Table II — Continued 



Element 


Source 


Temperature 
°K 


X 






1235 


0.0235 






1235 


0.0232 






1247 


0.0279 






1251 


0.0251 






1251 


0.0244 






1273 


0.0280 






1273 


0.0330 






1273 


0.0292 






1296 


0.0369 






1300 


0.0358 






1300 


0.0368 






1329 


0.0454 






1373 


0.0645 






1373 


0.0597 


Lead 


This work 


1323 


0.0022 






1349 


0.0025 






1373 


0.0036 






1423 


0.0046 






1474 


0.0076 






1475 


0.0072 






1523 


0.0112 


Arsenic 


Klemm & Pirscher (Ref. 38) 


1346 


0.60 






1483 


0.65 






1512 


0.70 


Antimony 


This work 


1073 


0.0151 






1073 


0.0143 






1170 


0.0239 






1174 


0.0243 






1272 


0.0461 






1371 


0.0835 






1473 


0.150 


Bismuth 


This work 


1273 


0.00262 






1373 


0.00544 






1423 


0.00832 






1473 


0.0113 






1548 


0.0193 


Zinc 


This work 


915 


0.00845 






925 


0.00962 






971 


0.0131 






973 


0.0168 






1017 


0.0208 






1020 


0.0200 






1024 


0.0247 






1069 


0.0321 






1071 


0.0325 






1072 


0.0394 






1117 


0.0486 






1119 


0.0526 
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Table II — Concluded 



Element 


Source 


Temperature, 
°K 


X 


Copper 


Hansen & Anderko (Ref. 41) 


1075 


0.300 




Rudolfi (Ref. 40) 


1103 


0.355 






1268 


0.424 






1336 


0.486 






1433 


0.598 






1503 


0.685 






1547 


0.763 


Silver 


Rudolfi (Ref. 40) 


1215 


0.168 






1307 


0.224 






1402 


0.300 






1513 


0.490 






1563 


0.622 






1592 


0.721 






1613 


0.794 






1623 


0.852 


Gold 


Hansen & Anderko (Ref. 41) 


643 


0.309 




di Capua (Ref. 43) 


1063 


0.439 






1303 


0.553 






1403 


0.637 






1448 


0.701 






1538 


0.775 






1568 


0.824 



3.8 Silicon-Tin, Silicon-Lead 

We have measured the solubility of silicon in tin, and our results are 
given in Table II and Fig. 9. The line corresponds to a = 8,145, b = 1.50. 
The measurements of Tamaru, 36 who used 92.5 per cent silicon (the 
principle impurities were iron and aluminum), give higher solubilities 
than we have found. 

Our measurements of the solubility of silicon in lead are also given 
in Table II, and the a-T plot (Fig. 9) shows the line, a = 19,830, b = 
4.58, drawn through the points. Moissan and Siemens 37 reported a num- 
ber of solubility measurements, all of which lie at appreciably lower 
silicon concentrations than those which we report. 

3.9 Silicon-Arsenic, Silicon- Antimony, Silicon-Bismuth 

Klemm and Pitscher 38 have obtained three points on the silicon-ar- 
senic liquidus curve that are listed in Table II. These three points 
have lead us to estimate the complete liquidus curve with the line, a = 
-49,990, b = -32.40, in the a-T plot of Fig. 10. 

Solubility measurements of silicon in antimony have been reported 
by Williams. 39 Our measurements, given in Table II, are in essential 
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Fig. 8 — a vs. T for silicon-aluminum, silicon-gallium and silicon-indium. 

agreement with these. However, we have used our measurements to 
establish the line that we believe to be a best approximation at present 
to the silicon-antimony liquidus curve. The line, a = 3290, b = —1.61, 
is shown with the experimental points in Fig. 10. 

Our measurements of the solubility of silicon in bismuth are recorded 
in Table II and plotted in Fig. 10. The line, a = 14,840, b = 2.06, has 
been drawn. 
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Fig. 9 — a vs. T for silicon-lead and silicon-tin. 
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Fig. 10 — a vs. T for silicon-arsenic, silicon-antimony and silicon-bismuth. 

3.io Silicon-Copper, Silicon- Silver, Silicon-Gold 

Rudolfi 40 has reported measurements of the solubility of silicon in 
copper. These are given in Table II, along with a value for the eutectic 
composition proposed by Hansen. 41 The a-T plot is shown in Fig. 1 1 
and the line, a = -11,910, b = -7.19, has been drawn. 

Arrivant 42 has measured the solubility of silicon in silver; his measure- 
ments are given in Table II. The line, a = —7910, b = —7.63, has been 
drawn through these data plotted as a versus T in Fig. 11. 

The solubility of silicon in gold has been measured by di Capua, 4 
and his measurements are listed in Table II. The plot of a versus T in 
Fig. 11 suggests that the data can be represented by the line, a = 
— 19,540, b = — 10.28. More weight has been attached to the four points 
at low temperatures than the high temperature points, where experi- 
mental error will lead to a greater scatter of points than at lower tem- 
peratures. 
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3.11 Silicon-Zinc, Silicon-Nickel 

Our measurements of the solubility of silicon in zinc, given in Table 
II, lead to an a-T plot, Fig. 12, which suggests the line, a = 4280, b = 
1.14. The measurements reported by Moissan and Siemens 37 lie at some- 
what lower values of the solubilities than those we report. 

The measurements of Iwase and Okamoto of the solubility of silicon 
in nickel, scaled from Hansen's plot, 1 are shown in Fig. 12 as an exam- 
ple of a binary silicon system for which the a versus T plot is clearly 
nonlinear. 



IV. DISCUSSION 

4.1 The Liquidus Curve Equation 

Equation (11) can be rearranged to give the following relationship: 

AH/ + a(l - x) 2 



T = 



ASi F - Rlnx + 6(1 - x) 2 ' 



(12) 



This will be called the liquidus equation since liquidus temperatures can 
be calculated as a function of x for each set of values of a and b. The val- 
ues of a and b used to approximate the available experimental data have 
been summarized in Table III. 
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Fig. 11 — a vs. T for silicon-silver, silicon-copper and silicon-gold. 
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Fig. 12 — a vs. T for silicon-nickel and silicon-zinc. 

The solid curves of Figs. 1 and 2 correspond to a = —150, b = f or 
the gemanium-gallium liquidus curve. The extrapolation of the a-T line 
to the melting point of germanium means that the liquidus curve ap- 
proaches the melting point of germanium with the limiting slope deter- 
mined by the heat of fusion and melting point of germanium. The ex- 
trapolation to lower temperatures of the a versus T line is in accord with 
the expected behavior of such liquidus curves — namely, that at low 
temperatures the log x versus 1/T plot will be linear. 

Table III — Coefficients a and 6 of the Liquidus 
Curve Equation for Various Binary Systems 



Element 


Germanium 


Silicon 


a 


6 


a 


6 




-5360 

-150 

1570 

(5700) 

1680 

8780 

-5600 

2640 

5505 

-7360 

-5500 

-4865 

4110 


-3.16 


0.56 

(1.90) 

1.08 

4.08 

-4.16 

1.98 

1.49 

-7.67 

-7.13 

-1.02 

1.75 


-4140 

3250 

11450 

(16600) 

8145 

19830 

-49990 

3290 

14840 

-11910 

-7910 

-19540 

4280 


-1.22 




0.83 




3.37 


Thallium 


(3.80) 


Tin 


1.50 




4.58 


Arsenic 

Antimony 

Bismuth 

Copper 


-32.40 

-1.61 

2.06 

-7.19 

-7.63 


Gold 


-10.28 




1.14 




— 
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The lilies in Figs. 3 through 11 have been drawn to indicate the com- 
plete liquidus curve and the amount of extrapolation of the experimental 
data required to give the full curve. The liquidus curves, plotted as T 
versus x or log x versus \/T, have not been included in this paper. These 
curves can be readily obtained from (12) and the constants of Table 
III. However, we have available complete sets of T versus x and log x 
versus \/T plots which will be supplied to anyone upon request. 

We have assumed that the distribution coefficients of the second com- 
ponents in germanium and silicon are small enough to have no influence 
on the liquidus curve. However, the solid solubility of gallium in germa- 
nium is high enough to cause us to expect that some influence should 
be expected. This is also true for aluminum-germanium. We will use the 
solid solubility measurement of Trumbore et al. 1 ' 45 to estimate the mag- 
nitude of the error introduced in the germanium-gallium system by 
ignoring the solid solubility. 

de Roche has reported that the liquidus curve points he has meas- 
ured that lie between 1.0 and 0.80 atom fraction of germanium lead 
to a heat of fusion of germanium of 8100 kilocalories. This calculation 
used ( 1 ) , and no correction was made for the solidus curve. If we accept 
the statement that the experimental liquidus curve follows ( 1 ) over the 
range 1.0 to 0.8 with A/// having the value 8,100, we can use the solid 
solubilities measured by Trumbore, et al. 1, and re-evaluate the heat of 
fusion of germanium. 

Equation ( 1 ) , modified to account for solid solubility, will be 

where x L is the atom fraction of germanium in the liquid phase and x s 
the atom fraction in the solid phase. The value of x L can be calculated 
using (1) with AH/ = 8100 and 1\° = 1210°K for various values of 
the temperature down to a value of x L of 0.8. The measurements of 
Trumbore et al. give x a at these temperatures. A plot of the logarithm 
of Xl/xs versus 1/T leads to a value of All/ of 7850 calories per gram 
atom. This value is essentially within the experimental error of 8100 ± 
200 cals estimated by de Roche. Consequently, the neglect of the solid 
solubility introduces a negligible error. 

4.2 Eutectic Temperatures and Compositions 

The eutectic temperatures and compositions that are consistent with 
the a and b values summarized in Table III are given in Table IV. The 
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Table IV — Eutectic Temperatures, T eu °C, Differences Between 
Melting Point of Second Component and Eutectic Tempera- 
ture (AT - T 3 ° - T eu ), and Eutectic Compositions in 
Atom Fractions of Germanium and Silicon 





Germanium 


Silicon 




r eu °c 


AT 


leu 


r eu °c 


AT 


*eu* 


Aluminum. . . . 

Gallium 

Indium 

Thallium . . . 
Tin 


424 
30 
157 
304 
232 
327 
736 
590 
271 
644 
651 
356 
398 
320 


236 
0.007 
0.2 
(0.3) 
0.9 
0.1 
78 
40 
0.05 
434 
310 
707 
22 
1.0 


0.28(0.30) 
5 X 10~ 6 
5 X 10-" 
(4 X lO" 4 ) 
0.003 
2 X 10"* 
0.58(0.59) 
0.18(0.17) 
2 X lO" 4 
0.36(0.36) 
0.25(0.26) 
0.27(0.27) 
0.044(0.055) 
0.002 


577 
30 
157 
304 
232 
327 
1073 
630 
271 
802 
830 
370 
420 


82 

6 X lO" 8 

1 x 10- 7 

(2 X lO" 8 ) 

4 X 10" B 

5 X lO" 8 
-259 

1.1 

3 X lO- 8 
281 
131 
693 

0.2 


0.120(0.121) 
5 X 10" 10 
2 X lO" 10 
(3 X 10-») 

1 x 10- 7 




9 X 10-'° 


Arsenicf 

Antimony .... 
Bismuth 

Silver 

Gold 


0.595(0.595) 
0.003 

1 x 10- 10 

0.32(0.30) 

0.125(0.154) 

0.31(0.31) 




4 X 10-* 


Cadmium 


— 



* Compositions in parentheses given by Hansen. 41 
t Not simple eutectic systems. 

eutectic temperatures have been taken from the literature and the 
corresponding composition calculated. There are a number of systems 
for which the eutectic temperatures are very near the melting point of 
the solvent. An estimate can be made of the eutectic temperatures and 
compositions by extrapolating the liquidus curve to the melting point of 
the solvent, then using this composition to calculate the freezing point de- 
pression, as was done in an earlier paper. 2 On the assumption that the 
solid solubility of germanium or silicon in the solvent element is negli- 
gible, the freezing point depression can be calculated from the following 
equation : 



A1 AHS eu ' 



(14) 



where .r eu is the atom fraction of germanium (or silicon) in the eutectic 
liquid phase, T 2 ° the melting point of the solvent, and AH 2 F its heat of 
fusion. Equation (14) follows from (1) applied to component 2. The 
freezing point depressions calculated in this manner, (see Ref. 14, Ta- 
ble B, pp. 284-305, for AH 2 F ), have been tabulated in Table IV, along 
with other eutectic compositions and temperatures. The germanium- 
gallium eutectic, calculated on the assumption of negligible solid solu- 
bility of germanium and gallium, is shown in the insert of Fig. 1. 
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4.3 The Relationship Between a and b 

In Fig. 13 we have plotted the values of b against a. It can be seen 
that, in general, the more positive the value of a, the more positive b 
becomes, and the more negative the value of a, the more negative b 
becomes. We have used this relationship between a and b to estimate 
the germanium-thallium and silicon-thallium liquidus curves. 

4.4 The Germanium-Thallium and Silicon-Thallium Liquidus Curves 

The available experimental measurements of the solubility of ger- 
manium in thallium are given in Table I, and the corresponding ct-T 
points are plotted in Fig. 3. The position of these points leads to a value 
of a of about 3500 calories per gram atom at 1075°C. On the assumption 
that the a-T relationship is linear and passes through this point, and that 
the value of a and b must be related in such a way as to fall on the curve 
of Fig. 13, we obtain for our estimate of the complete germanium- 
thallium liquidus curve, a = 5700, 6 = 1.90. 

From the position of the silicon-lead and silicon-bismuth a-T lines, 
we estimate a value for a of 12,500 at 1450°K for thallium. Then, pro- 
ceeding as for the germanium-thallium system, we conclude that a 
reasonable estimate of the a-T line for thallium would be a = 16,600, 
b = 3.80. The values of a and b for the thallium systems have been in- 
cluded in Table III in parentheses, and the eutectic compositions and 
temperatures are shown in Table IV. The a-T line for germanium- 
thallium is shown in Figure 3. 
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4.5 Implications of the a = a — bT Relationship 

In the preceding discussion we have used the a-T plots to estimate 
the complete liquidus curve from limited solubility measurements. We 
conclude that, within experimental error, the liquidus curves of a rather 
large number of binary systems for which germanium or silicon is the 
solid phase can be described by a two-parameter equation. Two systems, 
germanium-zinc and silicon-nickel, quite clearly fail to fall into the same 
class as the other liquidus curves. 

In Section 3.1 it was shown that, if the thermodynamic properties of 
the liquid phases are given by (9) and (10), it would follow that a 
would be a linear function of T along a liquidus curve. One cannot draw 
the conclusion, however, that, because a appears to be a linear function 
of T, the thermodynamic properties of the liquid phase must be given 
by (9) and (10). The partial molar heat and excess entropy are related 
to the partial molar excess free energy by partial derivatives at constant 
composition. The liquidus curves have given partial molar excess free 
energies along a line in A/Y, x, T, space. It is not possible to obtain 
partial derivatives of A/Y at constant composition from such limited 
information. The linear relationship between a and T only suggests that 
(9) and (10) describe the thermodynamic properties of these solutions. 
At present there is no information available which can be used as a 
direct check of this suggestion. 

From (3), (4) and (11) we can conclude that the following relation- 
ship is true along the liquidus curve if a is a linear function of T : 

AFx" = (a - 6T)(1 - a;) 2 . (15) 

If this relationship applied off the liquidus curve as well, it would follow 
from the Gibbs-Duhem equation that, for the other component, 

A/y = (a - bT)x 2 (16) 

and, for the integral molar excess free energy of mixing, 

AF m e = (a- bT)x(l - x). (17) 

It is possible for us to test the applicability of (16) to several of the 
systems studied by calculating the value of the activity coefficient of 
the second component at the eutectic temperature and composition. We 
can use the eutectic composition and temperature of a number of systems 
to obtain an experimental value of A/Y, which we can write in terms of 
the activity coefficient, t> , since 

A/Y = RT In 72 • (18) 
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Table V — Activity Coefficient of Component 2 at the Eutectic 



System 


^i 


7 , calculated 


Aluminum-Germanium 

Antimony -Germanium 

Gold-Germanium 


0.84 
1.08 
0.36 
0.99 
1.01 
0.41 


0.84 
1.02 
0.78 
0.97 




1.00 




0.38 







We have used (3), modified to account for a nonzero value of ACp (Cp 
of the liquid was assumed to be constant) as in (2), and applied it to 
the liquidus curve of the second component at the eutectic. The values 
of 72 have been obtained in this way for a number of systems, using the 
tabulated data of Kubashewski and Evans, 12 and are given in the first 
column of Table V. In the second column are values of 72 calculated from 
(16) and (18) using the eutectic temperatures and compositions in Table 
IV. 

A comparison of the activity coefficients in Table V reveals a fairly 
good agreement between the experimental activity coefficients and those 
obtained using the appropriate values of a and b in Table III.* An addi- 
tional bit of information is available for the silicon-silver system. Scha- 
del et al. have measured the vapor pressure of silver over silicon-silver 
melts at several temperatures in the range of 1000°C. The positive de- 
viations from ideality found by Schadel et al. are in fair agreement with 
the values we calculate using the a's and 6's found from the solubility of 
silicon in silver, but the temperature coefficient of the activity coeffi- 
cient they found is opposite to the one we predict. However, it is prob- 
able that not much weight can be attached to the experimental tempera- 
ture coefficients in view of the authors' comments on experimental 
errors. 

4.6 The Silicon-Thallium, Silicon-Lead and Silicon-Bismuth Binary 
Phase Diagrams 

We have used the a versus T linear relationship to fit the full liquidus 
curves of the silicon-lead and silicon-bismuth systems. When the values 
of a and b of Table III are put into (12), the liquidus curve equation, 
the curve rises from low temperatures to a maximum, passes through a 
point of inflection at x = 0.5 and then through a minimum before reach- 

* The germanium-silver s.ystem was not included in Table V because the solid 
solubility of germanium in silver at the eutectic is large and the activity coeffi- 
cient of silver in the solid solution is not known. 



GERMANIUM AND SILICON LIQUIDUS CURVES 



199 



ing the melting point of silicon, with the appropriate limiting slope. This 
behavior is in keeping with the fact that a miscibility gap is known to 
exist in these systems. If we assume that the values of a and b obtained 
from the liquidus curve at low temperatures give a reasonable estimate 
of the free energy of the liquid phases of the silicon-lead and silicon-bis- 
muth systems, and that our estimate of the a-T line for silicon-thallium 
is correct, we can calculate the complete phase diagram. 

The thermodynamic conditions for two liquid phases to be in equilib- 
rium lead to the following equation for the miscibility gap when the 
molar excess free energy of the liquid phases is given by (17): 

o(l - 2x') 

(19) 



T = 



R\n 



- 6(1 - 2x) 



The temperature at which this curve intersects the curve given by (12) 
is the so-called monotectic temperature, where two liquid phases, one 
rich in silicon and the other dilute in silicon, are in equilibrium with 
pure solid silicon. The maximum in the miscibility gap occurs at x = 
0.5. The temperature of the maximum in the miscibility gap, T c , the 
critical temperature, cannot be obtained directly from (19), but can be 
shown to be given by the equation 



T e = 



2R + b' 



(20) 



The critical temperatures predicted for the thallium, lead and bismuth 
systems are given in Table VI, along with the monotectic temperatures, 
T m , and monotectic compositions, x m ' and x m " . 

4.7 Liquidus Curves and Critical Temperatures 

The liquidus curves for three systems, silicon-indium, silicon-tin and 
silicon-antimony, are rather flat in the region of x = 0.5, when T is 



Table VI — Estimated Critical Temperatures, Monotectic Tem- 
peratures and Monotectic Compositions for Three Silicon 
Binary Systems with Liquid-Liquid Miscibility Gaps 



Element 


T c , °C 


T m , °C 


■V'm 


X"m 


Thallium 


1S62 
2047 
21S7 


1387 
1397 
1393 


0.0G 
0.03 
0.04 


0.94 


Le:ul 

Bismuth 


0.97 
0.96 
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plotted as a function of x. This indicates that the miscibility gap occurs 
at temperatures just below the hquidus curve in the supercooled liquids. 
If the free energies of the liquid phases were reasonably well approxi- 
mated by (17), the calculated critical temperatures would all lie below the 
liquidus curve temperature at x = 0.5. The comparison is made in Ta- 
ble VII for these three systems. Included in the table are the critical tem- 
peratures and liquidus temperatures for the three germanium systems, 
germanium-thallium, germanium-lead and germanium-bismuth, which 
have the highest calculated miscibility gaps. In every case the liquidus 
curve lies above the calculated critical temperature of the miscibility 
gap. This shows that the use of (17) for the free energy of the liquid 
phases is consistent with the presently known properties of these sys- 
tems. 

The systems germanium-arsenic, -copper, -silver and silicon-arsenic, 
-copper, -silver, -gold all have critical temperatures calculated from (20) 
that are higher than the temperature of the liquidus curve at x = 0.5. 
It can be shown, however, that this critical temperature is a lower one; 
that is, it is predicted that, above this temperature, the liquid phase will 
split into two liquids and below it there will be homogeneity. In all other 
systems, no critical temperature is predicted. 

Some indication of the degree of approximation to be expected from 
the above treatment of binary systems with miscibility gap can be ob- 
tained by applying this treatment to the copper-lead system, which has 
a miscibility gap that has been measured. We have evaluated a from 
the measurements of the solubility of copper in molten lead reported by 
Kleppa and Weil. A plot of a versus T can be approximated by the line, 
a = 6160, b = 0.85, although the experimental points clearly depart 
from this line near the monotectic temperature. Using these values of 
a and b in (12) and (21), the complete phase diagram was calculated. A 
comparison between certain critical compositions and temperatures ob- 



Table VII — Comparison of Liquidus Curve Temperatures 
and Calculated Critical Temperatures 



System 


T (liquidus, x = 0.5) 


T c (calculated) 




1114 
1132 
1122 

1590 
1581 
1583 


972 


Germanium-Lead 

Germanium-Bismuth 

Silicon-Indium 

Silicon-Tin 

Silicon-Antimony 


1090 
1007 

1563 
1487 
1213 
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Table VIII — Comparison of Calculated and Experimental 
Phase Diagram Features for the Copper-Lead System 



Critical temperature, T c ■■■ 

x e 

Monotectic temperature, T„ 

*»' 

x m " 



Calculated 



1263°K 
0.5 

1190 
0.26 
0.74 



Experimental 



1277°K 
0.35 

1227 
0.15 
0.67 



tabled by this calculation and those measured experimentally 41 can be 
made. This comparison is shown in Table VIII, where it can be seen that 
a fairly good approximation has resulted. 

4.8 The Thermodynamic Properties of Liquid Alloys 

Darken and Gurry 47 have used the a function [RT In 7/(1 — re) 2 ] to 
describe the thermodynamic properties of liquid alloys at constant tem- 
perature. The temperature dependence of a was not considered. Guggen- 
heim (lief. 15, pp. 83-84) has pointed out that the excess molar free 
energy of the carbon tetrachloride-benzene and carbon tetrachloride- 
cyclohexane systems over the temperature range from 30 to 70°C can be 
quite satisfactorily accounted for by an interchange energy, w, which is a 
linear function of temperature. This is equivalent to finding that a is a 
linear function of temperature. 

We have evaluated a along the liquidus curves of a number of other 
binary systems for which germanium or silicon is not one of the com- 
ponents. These systems have long liquidus curves and negligible solid 
solubility. The general features of the liquidus curves, on a T versus x 
plot, are the same as those of germanium and silicon binary systems. The 
a versus T plots of these liquidus curves* were nonlinear in every system. 

Scatchard 48 has discussed the relationship between the temperature 
dependence of a (in essence) and the volume change on mixing. The 
conclusion has been reached 49 that a volume change upon mixing that 
is greater than additivity indicates a positive excess entropy of solution, 
while a volume contraction indicates a negative excess entropy. Prigo- 
gine 50 and his co-workers have discussed this relationship in some detail 
in terms of their average potential model. (See also Ref. 51.) Exceptions 
can be expected to the above generalization, and such have been ob- 
served, but, in general, a close relationship between the excess volume 
and excess entropy of mixing can be expected. The correlation found here 

* Silver-tin, -lead; zinc-gallium, -indium, -tin, -cadmium. 
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between a and b (Fig. 13) suggests, then, that volume expansion upon 
mixing occurs for positive a and b while contraction occurs for negative 
values of a and 6. 

Kleppa 52 - 53 ,54 and Wittig 55 have measured the heats of mixing of quite 
a number of binary liquid alloys. Trends in heats of mixing with position 
of one of the components in the periodic system, keeping the other com- 
ponent fixed, were observed. Some of these trends can be satisfactorily 
accounted for" 2 in terms of Friedel's alloy solution model, which, among 
other things, takes into account the difference in valence expected for the 
two components. The relative valence effect does not seem to appear, 
in any simple way, at least, in the system we have studied. 

Wittig 86 has found some general trends in the heat of mixing that he 
has related to the position of the variable component in the periodic 
system. Germanium occupies a position as a variable component in the 
several series he has investigated. The heats of solution we have deduced 
from our analysis of liquidus curves do not appear to violate the trends 
Wittig has observed. 

V. SUMMARY AND CONCLUSIONS 

A two-constant equation, (12), has been used to fit the liquidus curves 
of a number of binary systems for which either germanium or silicon is 
one of the components. This equation gives a limiting slope of the liquidus 
curve at the melting point that is consistent with the known heats of 
fusion of germanium and silicon and the known solid solubilities. The 
form of the equation at low temperatures is also consistent with the ex- 
pected thermodynamic properties of the two phases. 

The general form of the equation suggests that the thermodynamic 
properties of the liquid alloys are simple functions of composition and 
temperature. It was pointed out that a first approximation to the rela- 
tive pai'tial molar heat of solution of germanium or silicon in the liquid 
alloys would be given by the equation 

AHi = a(l - x? (21) 

and a first approximation to the relative partial molar excess entropy 
by the equation 

A& e = 6(1 - x)\ (22) 

where a and b arc independent of temperature and composition. This 
means that the excess free energy would be given by 

A/Y - RT In yi = (a - bT)(l - x)\ (23) 

If solutions having the thermodynamic properties given by these three 
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equations were in equilibrium with pure solid germanium or silicon, the 
liquidus curve equation would have the form 

T = ***' ±-^L= x)2 (24) 

ASi" - Rkix+ 6(1 - x) 2 ' 

This is the same equation found empirically. 

The constants a and b of (24) were obtained from the straight line 
drawn through the points representing each solubility measurement 
when plotted as the function a versus T: 

= RTln 7l = TAS/ - AH/ - fi In *> = _ bT ( } 

(l - xy (1-.-P) 2 

It cannot be concluded however, that the thermodynamic properties 
of the liquid alloys must be given by (21), (22) and (23) as a result of 
the observed linearity of a with T. The linearity of a as a function of T 
shows that (23) gives the activity coefficients of germanium or silicon 
along the respective liquidus curves. It was found, however, that the 
activity coefficient of the other component at the eutectic in several 
binary systems was fairly well approximated by the equation resulting 
when it was assumed that (23) applied at all temperatures and composi- 
tions, and the Gibbs-Duhem relationship was used. Based on this as- 
sumption, the extent of the miscibility gaps in the silicon-thallium, -lead, 
-bismuth binary systems were estimated. 

The two-constant liquidus curve equation, (24) , can be used to provide 
a best present estimate of the complete liquidus curves for 26 different 
binary systems for which germanium or silicon is one of the components. 
The constants for each system are given in Table III. Sets of the liquidus 
curves, plotted as T versus x or log x versus 1/T, may be obtained from 
the authors upon request. 
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